Relativistic effects on information measures for hydrogen-like 

atoms 



K. D. Ser 

School of Chemistry, University of Hyderabad, Hyderabad 500046, India 

Jacob Katriejll 
Department of Chemistry, Technion, Haifa 32000, Israel 
(Dated: February 23, 2009) 

Abstract 

Position and momentum information measures are evaluated for the ground state of the rela- 
tivistic hydrogen-hke atoms. Consequences of the fact that the radial momentum operator is not 
self-adjoint are explicitly studied, exhibiting fundamental shortcomings of the conventional uncer- 
tainty measures in terms of the radial position and momentum variances. The Shannon and Renyi 
entropies, the Fisher information measure, as well as several related information measures, are 
considered as viable alternatives. Detailed results on the onset of relativistic effects for low nuclear 
charges, and on the extreme relativistic limit, are presented. The relativistic position density de- 
cays exponentially at large r, but is singular at the origin. Correspondingly, the momentum density 
decays as an inverse power of p. Both features yield divergent Renyi entropies away from a finite 
vicinity of the Shannon entropy. While the position space information measures can be evaluated 
analytically for both the nonrelativistic and the relativistic hydrogen atom, this is not the case for 
the relativistic momentum space. Some of the results allow interesting insight into the significance 
of recently evaluated Dirac-Fock vs. Hartree-Fock complexity measures for many-electron neutral 
atoms. 
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I. INTRODUCTION 

The celebrated Heisenberg uncertainty principle of quantum mechanics is specified 
by means of the position and momentum variances, that are defined in terms of the expec- 
tation values of the corresponding (hermitian) operators. Several authors have pointed out 
that for bimodal distributions this formulation does not provide an adequate measure of 
the uncertainty of the measurable involved. A variety of information measures have been 
proposed and investigated in a fairly broad range of contexts. The most familiar of these 
are due to Shannon [sl and Fisher 4]. They are being increasingly applied in studying the 
electronic structure and properties of atoms and molecules, and play an important role in the 
rapidly developing field of quantum information and its anticipated technological offspring. 

Due to its fundamental importance in natural sciences, the hydrogen atom has been 
extensivelystudied from the information theoretical view point. Consequently, the Shannon 
entropies |^, 6, 7| and Fisher information fsl of the non-relativistic hydrogen-like atom have 
been studied in considerable detail. We note, with some surprise, a glaring omission, the 
relativistic hydrogen atom, formulated in terms of the Dirac equation, which is the subject 
of the present study. Future studies of relativistic effects on the information measures in 
many-electron atoms will certainly benefit from the presently derived results. Indeed, very 
recent work on Dirac-Fock vs. Hartree-Fock complexity measures for neutral many-electron 
atoms allows certain comparisons to be made with results obtained in the present article, 
that shed additional light on the significance of that study. 

The Shannon information entropy Sr of the spatial electron density p(r) is defined as 



Sr = - J p(r) In p(r) dr , (1) 
and the corresponding momentum space entropy 5*^ is given by 

s, = - [ n(p)inn(p)rfp, (2) 



where n(p) denotes the momentum density. The densities p(r) and n(p) are each normal- 
ized to unity and all quantities are given in atomic units. These two densities are obtained 
from the corresponding position and momentum space wavefunctions, that are the Fourier 
transforms of one another. The Shannon entropy sum St = Sr + Sp contains the net infor- 
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mation and obeys the well known lower bound derived by Bialynicki-Birula and Mycielski 

St = Sr + Sp>n{l + \n7r) , (3) 

where n is the number of dimensions. The lower bound is attained by a Gaussian distribu- 
tion. This entropic uncertainty-like relation represents a stronger version of the Heisenberg 
uncertainty principle of quantum mechanics. The individual entropies Sr and Sp depend on 
the units used to measure r and p respectively, but their sum St does not, i.e., it is invariant 
under uniform scaling of coordinates. 

The Shannon entropies provide a global measure of information about the probability 
distribution in the respective spaces. A more localized distribution yields a smaller value of 
the corresponding information entropy. For applications of Shannon information entropy in 
chemical physics we refer the reader to the published literature 

In the context of the quantum theory of one-particle systems the Fisher position infor- 
mation measure is defined as 

-/^- 

and the corresponding momentum space measure is given by 

Ip = j^^dp. (5) 

For a general definition of the Fisher information measure and a careful exposition of its 
significance we refer to the definitive monograph by Rao [12]. 



The individual Fisher measures are bounded through the Cramer-Rao inequality 1^, 
13j according to > — and Ip> —, where V^'s denote the corresponding spatial and 

Vp 

momentum variances, respectively. In position space the Fisher information measures the 
sharpness of the probability density, and for a Gaussian distribution is exactly equal to 
the inverse of the variance [l4|. A sharp (smooth) and strongly localized (well spread-out) 
probability density gives rise to a larger {smaller) value of the Fisher information in the 
position space. With a differential probability density as its content, the Fisher measure 
is better suited to study the localization characteristics of the probability distribution than 



the Shannon information entropy 15|, [iGj. Unlike Sr + Sp, for which eq. [3] specifies a lower 



bound, general bounds are as yet unknown for the Fisher product Ijlp. Since localization 
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(i.e., low uncertainty) means high values of the Fisher information measures, the counterpart 
of the Heisenberg or Shannon bound should be an upper bound on the product of position 
and momentum Fisher information nieasures. For a single particle under the influence of 
a central potential Dehesa et al. |J/7(] have very recently reported a lower bound on the 



Heisenberg product 
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2l| which can be directly related to the Fisher information. 



For the application of the Fisher information measure as an underlying guideline for the 
formulation of fundamental physical principles we refer to the recent book by Frieden fl^ . 
and for applications to the electronic structure of atoms, to the pioneering work of Dehesa 



et al. 
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25|. 



A widely used generalization of the Shannon entropy is the Renyi entropy. The Renyi 
position entropy (that, when a more precise designation is required, we shall address as the 
a- Renyi position entropy) is defined as [26 1 



1-a 



log / [p{r)Y4:'Kr dr 



The symmetrized Renyi position entropy is 



(6) 



(7) 



where 



1 -s- 



1 + 



1 , 11 

-1 < s < 1, i.e., - < a, < oo, — 1-7 = 2. 

2 a b 



The Renyi momentum entropies are similarly defined in terms of the momentum density 
n(p). For a = 1 one obtains, using I'Hopital's rule. 



lim Hi'^ 



p{r) log[p(r)]47rr^(ir = Sr 



where Sr is the Shannon position entropy [3|. For an n dimensional system, the sum of 
the a-Renyi position entropy and the 6-Renyi momentum entropy, for a and b satisfying 
^ + 1 = 2, was recently shown by Bialynicki-Birula ^] to satisfy the inequality (uncertainty- 
like relation) 



n 



(8) 
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The properties of the Shannon entropies and several other information measures under 



coordinate scaling have been examined in ref. 



27| . It was pointed out that upon scaling the 



coordinates via r = (r, normalization of the density requires that 



It follows that 



1 fr 



= (C'^'"''^ [pix)rAnx'dx^ = 3 log(C) + 



(r) 
a 



The scaling of the coordinates introduced above implies scaling of the momenta according 
to p = 2, so 

i7lf) = -31og(C) + i/5'^ 

Hence, 

Note that a and a' are entirely independent of one another. 

For a system whose hamiltonian is of the form Ti = T + \V{r) , where T = — ^V^, if 
the potential is homogeneous, i.e., V{(^r) = C^y(r), scaling of the coordinates via f = is 



equivalent to scaling the coupling constant via A = C^"'"^A [271]. Hence, for such potentials 
the sum of an a-Renyi position entropy and an a'-Renyi momentum entropy is independent 
of the coupling constant A. 

Since the position density has dimensions of inverse volume and the momentum density 
has dimensions of inverse momentum cubed, the Shannon, Fisher and Renyi entropies can 
be converted into quantities that have dimensions of length or momentum. This property is 
used to facilitate comparison among the different information measures. We refer to these 
transformed quantities as "length" and "impetus" , respectively. This practice allows a clear 
distinction between variance-based uncertainty measures and information based measures 
that have units of position or momentum although they do not involve expectation values of 
the corresponding quantum-mechanical operators. "Impetus" is a pre-Newtonian synonym 
of momentum. 
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A different class of information measures, involving ratios between relativistic and nonrel- 



ativistic densities 



, |29| , is considered as well. Since these information measures are pure 
numbers they cannot be transformed into quantities with units of length or momentum. 
However, they allow the onset of relativistic effects upon increase of the nuclear charge to 
be followed very transparently. 



A third class of information measures is represented by the Tsallis entropy SOj, 

/>oo 



a — I 

that has been invoked for non-extensive systems, and which is not homogeneous under 
coordinate scaling. The Tsallis entropy is closely related to the Renyi entropy via 

e, = -l^(l-exp((l-a)i7W)) . 

For a —>■ 1 these two entropies coincide with one another as well as with the Shannon entropy. 

Upon attempting to evaluate the various uncertainty and information measures for the 
relativistic hydrogen atom we encountered three somewhat surprising obstacles. The first 
has to do with the fact that the radial momentum does not have a proper (self-adjoint) 
quantum- mechanical counterpart. This fact has been known for a while, and its conse- 
quences in the present context are explained below. Another surprise had to do with the 



fact that t 
versy 



31 



le momentum-space solutions of the Dirac equations are still subject to contro- 



32l |. Finally, we find it intriguing that for the relativistic hydrogen atom the 
position space information measures could be evaluated analytically, but for the momentum 
space measures we had to apply numerical integration. The first difficulty suggests that the 
various uncertainty-like principles that do not involve the variance of position or momentum 
have even stronger merits for multi-dimensional systems than those extensively pointed out 
in previous studies on one dimensional systems. We addressed the second difficulty by ad- 
hering to the version of the momentum wavefunction that in our judgement was the most 



natural and straightforward [33| without bothering to explore its equivalence (or lack of it) 
with other formulations. 

This article is structured as follows: In sections 2 and 3 we consider the radial position and 
momentum variances for the nonrelativistic and the relativistic hydrogen atom, respectively, 
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presenting very explicitly the consequences of the non self-adjointness of the commonly 
invoked radial momentum operator. The Shannon entropies are investigated in section 4, 
and the Fisher information measures in section 5. In section 6 we study the Renyi entropies, 
also considering the average density, that is essentially a special case of these. This allows a 
brief discussion of complexity measures, with a rather surprising comparison with a recent 
Dirac-Fock study of neutral many-electron atoms jsl . Scale invariant entropies are discussed 
in section 8, and some concluding remarks are made in section 9. 

II. UNCERTAINTY MEASURES FOR THE NONRELATIVISTIC 
HYDROGEN-LIKE ATOMS 

The most familiar measures of uncertainty are the position and momentum variances, 
which for a one dimensional system defined over the whole real axis are given by 
{5x) = < x'^ > — < X and {6p) = a/< > — < p where 



Here, iplx) and (j){p) are the position and momentum wavefunctions, which are the Fourier 
transforms of one another. Eq. [9] emphasizes the correspondence between position and 
momentum space expectation values of the (hermitian) position and momentum operators, a 
correspondence that, as we shall explicitly demonstrate below, fails for the (non self-adjoint) 
radial momentum operator. 

A. Position and momentum variances for spherically symmetric three-dimensional 
systems 

Qiang and Dong [34], following a time- honored tradition, suggest that the radial momen- 
tum operator in the coordinate representation, using atomic units in which h = 1, is 
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They justify this expression by noting that 



2 - _ / ^ 2 9_ 



is (up to a multiphcative constant) the radial part of the Laplacian (i.e., the kinetic energy 
operator). 

Note, however, that the n dimensional generalization 

d n — 1 



satisfies 



dr 2r 



. n-l d in - l)(n-3) 
-p = \ h -^^ 



which, for n = 2 and n > 3 does not agree with the radial part of the Laplacian, 



Straightforward integration by parts yields 



'o 



d n — 1 
dr 2r 



r'^-^dr = 



(10) 



n 

where Sn = ^7vy is the surface area of the n-dimensional unit sphere. 



Paz 



35| has recently rigorously shown that pr is not self-adjoint, and has no self-adjoint 



extension (c/., also, [36|, l37|, l38|, for earlier discussions of this issue). This fact, which is 
more closely considered below, suggests that using the variance of the radial momentum 
as an uncertainty measure for (spherically symmetric) three-dimensional systems may be 
questionable. 

B. Position space wavefunction of the nonrelativistic hydrogen atom 



The ground state nonrelativistic hydrogenic wavefunction V'o,o,o(^) 
yields the density 



PNR{r) 



TT 



exp(— 2Zr) . 



' exp(— zrj 
(11) 
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The expectation values < ^ >= ^ and < >= -p- yield 

. , 2 2 ^ 0-8660 

[or) = v< > ~ < r >^ = — ^ — - — . 

Similarly, < "Pr >= and < >= Z^. The latter value is equal to < — >, because, for 
a spherically symmetric wavefunction, the angular part of the Laplacian makes a vanishing 
contribution. 

Using these position space expectation values to evaluate (5pr) = "Z, we obtain the 
uncertainty product {br){bp,^ ~ 0.8660. 



C. Momentum space wavefunction of the nonrelativistic hydrogen atom 

The momentum space wavefunction for the ground state of the nonrelativistic hydrogen 
atom [391] 



Xo,o,o(p) 



1 /2 



V 



+ 1 



yields the momentum density 



(12) 



-\ -4 



(13) 



in terms of which we obtain 



<p>= I pUNR{p)Anp'^dp = — ^ 0.848826Z 
'0 Svr 



and 



< P^ >= I p IlNR{p)4:iTp^dp = Z^ 

'0 



The latter value agrees with the position space expectation value of — V^, but the former 
does not agree with the position space result, eq. [10 

Using these momentum space expectation values we obtain 



Z 



(Sp) = — V97r2 _ 64 ^ 0.5287Z . 
on 
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Along with the value of 5r obtained above we get 



(5r)((5p) ^ 0.4578 , 

which is less than |. This is probably a manifestation of the questionable status of the radial 
momentum, pointed out above. 

Messiah j40(] shows that, in one dimension, < >< p^. > > ^. For a spherically 
symmetric system < >= 3 < > and < >= 3 < p^ >. Hence, < >< > > |. 
The results quoted above imply that for the nonrelativistic hydrogen atom < >< p'^ >= 3, 
which is larger than the lower bound derived by Messiah. 

We conclude this section by emphasizing that the operators and p"^ can be expressed 
(in Cartesian coordinates) in terms of manifestly self-adjoint operators. This is not the case 
for r and p. Further consequences of this distinction are presented in the following section. 



III. UNCERTAINTY MEASURES FOR THE RELATIVISTIC HYDROGEN-LIKE 
ATOMS 

A. Relativistic position uncertainty 



For the ground state of the Dirac hydrogenic atom the (spin up) wavefunction is of the 
form 



= Nr^-^ exp{-Zr) 



(g \ 



\wYiJ 



(14) 



where G = a/1 + 7, g = y/1 — 7, Yq = cos{6) and Yi = sin(6') exp(z(/ 
factor is given by = (2^) where 7 = [1 — (ZaY]^. Here a 

^ ,/87rr(27+l) ' ^ V 7 J 



it follows that 7 — > 1 and we obtain the nonrelativistic wavefunction y — exp(— Zr). 
The ground state position wavefunction yields the position density 



The normalization 

137.03600 

structure constant. In the limit Z ^ or c ^ 00 [remembering that in atomic units a = 1/c] 



27+1 



47rr(27 + 1 



-r 



2(7-1) 



exp(— 2Zr) 



(15) 
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In the extreme relativistic limit Z - the position density obtains the form 41 1 



pER[r) = - — r exp 

zna \ a 

Using the relativistic wavefunction eq. (fT4|) . we obtain 

27 + 1 

and 



hence, 



{6r) 



2Z 



For future reference we define the ratio P = \J ^5^^, which is plotted in Fig. 1. Evaluating 
< Vr >R= and 

< pI >r= ~< pI >nr (^{aZf + h^aZY + ■ ■ ■ ^ (16) 

it follows that [b^r) = '^/^Tf- Hence, 



A singularity is observed for 27 — 1 = 0, that yields Z = ^ ^ 118.68. We do not know 
what significance to assign to this nuclear charge. 

Since the small components of the Dirac wavefunction for the hydrogen atom depend on 
the angular coordinates, the expectation value of the Laplacian is not the same as that of 
Cr, in spite of the fact that the ground state density is spherically symmetric. The angular 
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part of the Laplacian yields 



;Y2^2{7-i)gxp(-2Zr)rfr MG, 0, -tgYo, -tgYi)\L'^ 



(G \ 




1-7 



27(27 - 1) 



Adding the value of < '^B\Vr\^D >, eq. [lEl we obtain 



7(27-1)- 



The singularity at Z ^ 118.68 remains 



B. Relativistic momentum uncertainty 



(17) 



We use the expression for the relativistic ground state momentum wavefunction due to 
Sheth 331] . Denoting the radial momentum variable by p and defining 



x = (7 + 1) arctani 



P 



(18) 



the momentum density can be written in the form 



r(7 + i: 



F{p) 



2Z'^TT2r{'y + i) (E)2 



:f)2 + i 



7+1 



where 



F{p) = (7 + 1) sm\x) + ^-^R{pY 

T 



(19) 



and R{p) = (7 + 1) cos(x) — ^ sin(x). T is the familiar F-function. 

In the nonrelativistic limit 7^1 this expression reduces to Il]\fji{p), cf. eq. f|T3|) . 

In the extreme relativistic limit, Z ^ the momentum density is obtained by using 
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L'Hopital's rule to evaluate lim^^o "f = 



. It is found to be of the form 



7=0 



27r2 [ap{l + {apY 



{apf + 1 



(ap H ) arctanfap) 

ap 



(20) 



The long range decay of the momentum density can be established by noting that for 
p — > oo 



1 



(f)'((f)' + l 



7+1 



P J 



(27+4) 



Furthermore, 



so that 



TT 



X ^ (7 + 1)2 



F ^ = (7 + 1) ( 1 + cos^ ( (7 + 1) 



TT 



r V -2, 

To evaluate < > we integrate over p numerically between p = and p = Pm, where Pm is 
chosen large enough for the density to be close enough to its asymptotic form, and add an 
integral over the asymptotic momentum density between p = pm and p = 00. Hence, 



< p" 



Jo v^r(7 + 2) 



-(27-1) 



27 - 1 



The ratio fi 



<P^>NR 



is plotted in Fig. 1. The expectation value < p^ >r was evaluated 



for the values of Z considered by Qiang and Dong [S^] • The results agree with those evaluated 
analytically by using the position space expectation value of the Laplacian, eq. (JT71) . for 
Z = 1, 11, 37, to ten decimal places, and for Z = 87 to eight decimal places. The values of 



y <p1 >, evaluated in terms of the ground-state wavefunction in the position representation, 
using Eq. f|T6|) , agree with the values of (Apr)_R in {s^]. It would be nice to show analytically 
that < > is indeed equal to the right hand side of Eq. f|T7|) . 

The results presented above clearly expose the difference between the self-adjoint Lapla- 
cian and the non-self adjoint radial momentum. 
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IV. SHANNON ENTROPIES FOR THE HYDROGEN ATOM 



A. Nonrelativistic position and momentum entropies 



The nonrelativistic position entropy S^^ — — 4Trr'^pNR{r) log (^PNR{r)Jdr can be eas- 
ily evaluated for the hydrogenic nonrelativistic ground state position density, yielding 



S^^ = 3 + log(7r) - 3 log(Z) ^ 4.1447299 - 3 log(Z) . (21) 



Similarly, the nonrelativistic momentum entropy S^^ = ATrp'^IlNji{p) \og \^NRip)j dp 
can be evaluated for the hydrogenic nonrelativistic ground state momentum density, yielding 

+ 5 log(2) + 2 log(7r) + 3 log(Z) ^ 2.4218623 + 3 log(Z) . (22) 



It follows that 



S^^ + S^^ = — + 5 log(2) + 3 log(7r) fa 6.5665922 
3 

> 3(l + log(7r)) ^ 6.4341897 . 



B. Shannon length and impetus 

Since exp(5'r) has dimensions of volume, we define the Shannon length Rs via 

47r 

—Rl = eMSr) . 

We similarly define the Shannon impetus Ps via 

47r 

—Pi = eMSp) . 
For the nonrelativistic hydrogen atom we obtain 



R 



1 

NR 



3\3 e 2.46972 

— I 

A Z Z 
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and 



hence, 



P/^ = 2(37r)^exp(-^)Z^ 1.39071Z, 

9 



i?^^P^^ ^ 3.43466 



From the entropic uncertainty-hke relation 

Sr + Sp>3{l+\n{7T)) 

it follows that 



3 



2 



RsPs > J ' vre ^ 3.28639 . 

This lower bound is, indeed, lower than the product of the hydrogenic Shannon length and 
impetus, but fairly close to it. 



C. Relativistic position entropy 

The position entropy can be evaluated analytically in terms of the relativistic ground 
state position density, eq. f|T5|) . to yield 

= log (^^^^f^) + (27 + 1) - 2(7 - 1)^(27 + 1) (23) 



^iog(r(^)) 



dz 



where "^{z) is the Digamma function, defined as "^{z) 

In the nonrelativistic limit 7 = 1 so r(27 + 1) = r(3) = 2 and reduces to S^^, eq. 

EU 

The relativistic correction is 

Si'-Sr = 2(7-l)(l-^(27 + l)) + log [^^^^) - -("^)'+ Q - T^) ("^)' + - ■ ■ ■ 

(24) 

In the extreme relativistic limit 

S-^ = log — - - 2^ + 1 ^ -14.463580 (25) 



16 



where C = 0.5772156649.. is Euler's constant (more commonly denoted 7, a notation we 
avoid for an obvious reason). 



D. Relativistic momentum entropy 



The relativistic momentum entropy was evaluated by means of numerical integration. To 



examine the relativistic corrections to the position and momentum entropies more closely, 

Eq. El yields ^ -1 - O.U75{aZf + ■■■ 



we consider cr^ = "(^^[2 and ap = -^j^ 



and a numerical fit yields dp ^ 1.80 + 0.65(aZ)2 -\ , so + dp ^ 0.80 + 0.2(aZ)2 + 



E. Relativistic Shannon length and impetus 



For the relativistic Shannon length we obtain 



R 



R 



— (3r(27+l)j exp 



l + ^(7-l)(l-M/(27 + i; 



Hence, 



[is 



R§ 
R^^ 



r(27 



exp 



^(7-l)(l-*(27 + l)) 



The ratio of the relativistic and nonrelativistic Shannon impeti can be obtained in terms of 
the numerically determined ap, as follows. 



pR 



{aZf 



-CTr, 



Hence, 



Psfis ~ 1 + 0.27(aZ)2 + 



(26) 



In Fig. 1 we present the ratios of the relativistic to nonrelativistic Shannon lengths and 
impeti, along with the corresponding ratios of the root mean square radius and momentum, 
and that of the Fisher lengths and impeti (to be discussed below). We note that for large Z 
the relativistic effect on the momentum uncertainty measures is larger than on their position 
counterparts. This is most pronounced for the root mean square position and momentum. 
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whereas the relativistic effects on the Fisher measures are almost symmetrical. This is 
most likely due to the fact that the Fisher measures are sensitive to the local oscillations 
of the distribution rather than to its long range behavior, where the relativistic momentum 
distribution varies the most. 



V. FISHER INFORMATION MEASURES FOR THE HYDROGEN ATOM 

The Fisher position information measure, Ir — /o°°47rr^^ (^)^^^ ^® easily evalu- 
ated for the ground state of the relativistic hydrogen-like atom, yielding 

= ^ =^<pi>^ ir (i + (o^zr + lio^zr + f +...), 

where I^^ = 4Z^. Since has a singularity when 27 — 1 = 0, i. e., aX Z k, 118.68, we do 
not examine the extreme relativistic limit. 

The Fisher momentum information measure Ip — /o°°47rp^jY^ (^) nonrela- 
tivistic momentum distribution can be evaluated analytically, yielding I^^ — For the 
relativistic momentum density the Fisher information was evaluated numerically. 



A. Fisher length and impetus 

The Fisher position information measure has dimensions of inverse area. We define 

_i _i 
the Fisher length as Rp = Ir ^ ■ Similarly, we define the Fisher impetus Pp = Ip . For the 

nonrelativistic hydrogen atom we obtain Rp^ — and Pp^ — The ratio between the 

relativistic and the nonrelativistic Fisher lengths is 

pi? 1 1 

(3p^-^^{2^-lY^l- -{aZf - -{aZf + • • • , 

The corresponding ratio of Fisher impeti is evaluated by fitting the numerically evaluated 

ratios to obtain 

pR 

liF = « 1 + 0.4166(aZ)2 + 0.23(aZ)^ + • • • , 

hence, 

pFl^F ~ 1 - 0.0834:{aZf + ■■■ . (27) 
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We note that the leading relativistic term is negative, unhke the corresponding term for the 
product of Shannon length and impetus, eq. [26l 

The ratios of the relativistic to nonrelativistic Fisher lengths and impeti are presented in 
Fig. 1. 



VI. RENYI ENTROPIES FOR THE HYDROGEN ATOM 



A. Renyi position entropies 



The Renyi position entropy is defined in eq. [61 The hydrogenic nonrelativistic ground 
state position density yields 



(28) 



For a — > 1 this expression reduces to h[^'^^^ = log (-p-) + 3, which is the well-known 
hydrogenic Shannon position entropy. Substituting in eq. [7] we obtain 



^(.-)=logf^UAlogfii±^ 



The hydrogenic relativistic ground state position density yields 

Hi-^^ = log(^^^||^)+2(7-l)log(a) 

'rf2(7- l)a + 3' 



(29) 



+ 



1 - a 



log 



r(27 + i) 



(27 + l)log(a) 



This expression is finite when the argument of the F-function satisfies 2(7 — l)a + 3 > 0. 
This condition can easily be traced back to the singularity of the relativistic density at the 
origin, cf. eq. [151 Since < 7 < 1, this condition holds for all Z when a < |. For a > | 
divergence will take place when Z > ^^^^^ ■ For 7 = 1 eq. [291 reduces to eq. [281 In the 
limit a — > 1 we obtain the relativistic Shannon entropy, eq. [231 For 7 — > this expression 
yields S^^. 
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The relativistic correction AHjf"^ = Ha'^^ — H^'^^^ can be expanded in the form 

2 )^ ' a-\ a-l^V r(27 + l) ) 

^^((aZr + (aZ)V4 + ---)+aQ-^) {o^Zf 



a log(a) 



Using this expansion we obtain 



In the extreme relativistic hmit 



and 



(30) 



Eq. [30] can also be obtained from eq. [29l by taking the limit 7 — > 0. For a — > 1 this 
expression yields eq. [251 

The results presented above imply the commutativity of the diagram 

PNR[r) =^ Ha > 

T T T 

pnir) =^ H^:'""^ 

where 



t stands for lim 
I stands for lim 

7-^0 

stands for lim 



1 /"°° 

X(r) =^ Y stands for Y = / 47rr^[X{r)]''dr , 

i.e., the fact that whenever more than one path (respecting the directions of the various 
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arrows) is available between any two nodes, the results along the different paths are identical. 



B. Renyi length 

Noting that exp{Ha^) has dimensions of volume we define the Renyi length Ra via the 
relation 

4.7T 

-Rl = exp(if(^)) . 



It follows that 



and 



Similarly, 



e 



limi?^«= - ^ = R 



1 1 _ 

) 



yielding -R^"^ in the limit 7 — >^ 1. 



C. Renyi momentum entropies 

The hydrogenic nonrelativistic ground state Renyi momentum entropy is 



r{p,NR) - / 7r^Z^\ 1 



Hr"" = log [-^) + YZTb (A^)) ' (31) 
where, 

m)-'-^fiy^.i)"V^y-^%fi. (32) 

The dependence of Ha'^^\ eq. EHl and H^'^^\ eq. [3ll on the nuclear charge Z is such 
that the sum, for any choice of a and &, is independent of as demonstrated above for 
arbitrary homogeneous potentials. 

Noting that J(l) = 1 we obtain, for 6-^1, the nonrelativistic Shannon momentum 
entropy, eq. [22l 

The nonrelativistic momentum density behaves, at large p, as IlNRip) ~ ^1 so that the 
integral in H^f' diverges unless 86 — 2>1, or6>|. Indeed, for 6 = | the numerator of 
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eq. [32] vanishes. Note, however, that this value of b is below the lower bound b > ^ allowing 
the definition of the symmetrized Renyi entropy, eq. [71 

The extreme relativistic momentum density, eq. [201 behaves, for p oo, like 

so that the integral in the expression for H^'^^^ behaves like ^^_2 ■ Hence, the integral 
diverges unless 46 — 2>lor6>|. 

The behavior of the relativistic momentum density is more subtle. For p ^ oo the variable 
X, defined in eq. [T8|, satisfies x — > (7 + As long as 7 < 1 one finds that F{p), defined 
by eq. [19], becomes a (7 dependent) constant, so that for large p the relativistic momentum 
density decays as n/j(p) ~ 2-^_i_4 . It follows tliBit tliG iiitcgrRiid in tliG GxprGSsioii for h'^'^'^ 
converges provided that (27 + 4)6 — 2>lor6> 3^^. For 7 = this expression yields 
6 > |, in agreement with the result obtained above for the extreme relativistic momentum 
density, but for 7 = 1 this expression yields 6 > |, which is larger than the bound 6 > | 
obtained above for the nonrelativistic momentum density. This is a consequence of the fact 

2 

that by taking the limit 7 — > 1 before the limit p ^ 00 one obtains F = 8^^^^, that for 
large p yields F ^ ^ rather than the constant obtained when the limits over 7 and p are 
taken in the opposite order. Since < 7 < 1, it follows that for 6 > | the Renyi entropy 
converges for all Z, for 6 < | it diverges for all Z, and for | < 6 < | it converges for 
Z < ^1 /(2 — — |- Substituting b = we obtain Z < ^^^^^ . Comparing with 

V a 

the results obtained above for h!^'^^ we conclude that h!^'^^ and h'^'^^ converge over the 
same range of Z when a and b are related via - + \ = 2. 

The relativistic Renyi momentum entropies can only be obtained numerically. The main 
point to note is that the various sums of Renyi position and momentum entropies exhibit a 
dependence on Z, unlike the nonrelativistic case. In Fig. 3 we show the sum of the Renyi 
position and momentum entropies, H^i^ + where a = and b = vs. s. The 
lowest curve corresponds to the lower bound presented in eq. [8], and the curve just above 
it is the nonrelativistic entropy sum. The relativistic entropy sums are all higher than the 
nonrelativistic one, exhibiting a rapid increase for higher s values. This behavior anticipates 
the approaching singularity of the relativistic Renyi momentum entropy for an appropriate 
value of s > |, that decreases with increasing Z, as clearly displayed in Fig. 3. Thus, for 
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Z = 100 the Renyi momentum entropy becomes singular for s = ^-^^ ~ 0.789. 

D. Renyi impetus 

The relation 

^P,' = exp(iff ) 

defines the Renyi impetus Pa. It follows that 

1 



pNR^l'^yZf 8 r(46 



4 J 2 Vv^ r(46) 

In the limit 6^1 this expression yields the nonrelativistic Shannon impetus. 

The relativistic Renyi impeti can be obtained from the numerically evaluated Renyi 
momentum entropies. 

From the uncertainty-like relation for the Renyi entropies, eq. M, it follows that the 
length-impetus product satisfies 

RaPb>ij^] a^b^ . (33) 

where ^ + 7- = 2. 

a b 

The ratios of the relativistic to nonrelativistic Renyi lengths {ISr} and impeti {^r) are 
presented in Fig. 2 for conjugate pairs (a, 6) = {(§, §)'(!' D' (|' D' d' D' (§' D' (§' §)}• 
Like the Shannon measures presented in Fig. 1, the Renyi impeti show a more pronounced 
relativistic effect than the corresponding lengths. The relativistic effect on the Renyi im- 
peti increases with increasing 6; the relativistic effect on the corresponding lengths (that 
correspond to decreasing a, satisfying ^ + ^ = 2) also increases, but more moderately. 

E. Average position and momentum densities 

The average position density is defined as 



< p >= J 47rr^^p(r)j dr , 
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While closely related to the 2-Renyi entropy, i.e., < p >= exp ^— ifg'^-' j , the average position 
density merits special attention since it has recently been invoked as a factor in a proposed 
measure of complexity 42, |43|. The average density is also known as the Onicescu informa- 
tion measure 4J], and is closely related to the linear entropy = 1— < p >, which is the 
q = 2 case of the Tsallis entropy [s^. 
For the hydrogen atom 

< Pnr >= ^ , 

and 

Z^TU-f - 1) 

<PR>^ 



7r247-2(r(27 + 1))2 ■ 

The relativistic expression becomes singular when 47 — 1 = 0, z.e. Z = ^ 132.68. The 
onset of relativistic effects is given by <pr>-<P'^R> ^ 0.055159(aZ)2 + 0.067737(aZ)^ + 
0.079452(aZ)6 + . ... 

The momentum density expectation value < 11 >= 47rj9^ ^n(p)^ dp for the nonrel- 
ativistic density was evaluated analytically, yielding < Unr >= ^g^a^a ~ o-208975 ^ rjj^^ 
relativistic counterpart, < >, can only be evaluated numerically. 

The leading terms in the Taylor series expansion of Z^{< II/j > — < Unr >) were 
obtained by differentiating the integrand, Z^Anp"^ (^(IIr{p)Y — (IliyR{p)Y^ , with respect to 
aZ, an appropriate number of times, evaluating it at Z = 0, and integrating over p. In this 
way we obtain 

Z'\< Ur> - < Unr >) ^ -0.254464(aZ)2 + 0.054446(aZ)^ + 0.001883(aZ)^ + ■ ■ • . 



or 



^^^^ ~ 1 - 1.21768(aZ)2 + 0.26054(aZ)^ + 0.00901(aZ)'^ + . 



< ^NR > 
F. Average length and impetus 

Noting that < p > has dimensions of inverse volume we define the average length Ra via 

—R\ =< p . 
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Similarly, the average impetus Pa is defined via 



As a consequence of the connection between the average entropy and the 2-Renyi entropy 
H2^^ the lengths and impeti related to these two information measures coincide. For the 
nonrelativistic hydrogen atom we obtain 

6i 1.81712 



and 

Furthermore 

Ha 



^ Z Z 



« = Z f'l^V ^ 1.04538Z 



11 / 



and 



hence, 



pNR 



1 - 0A6210{aZy - 0.14040(aZ)^ - 0.07719(aZ)^ + ■ 



1 + 0.40589(aZ)2 + 0.24265(aZ)^ + 0.16806(aZ)^ + 



1 - 0.0562(aZ)^ 



This expression is rather similar to the corresponding ratio for the Fisher length and impetus, 
eq. m 



G. Complexity measures 



The results presented above concerning the relativistic effects on the various information 
measures for the H-like atoms allow the evaluation of the statistical measure of complexity 
C, defined by Lopez- Ruiz, Mancini, Calbet (LMC) 42|, |43|]. The LMC measure C is given 
by 

C = H -D , (34) 
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where H denotes a measure of information and D represents the so called disequilibrium or 
the distance from equilibrium (most probable state). The form of C is designed such that it 



vanishes for the two extreme probability distributions corresponding to per 
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ect order (H = 0) 
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47|, 



48|, that 



and maximum disorder (D = 0), respectively. It is only very recently 
the studies on the electronic structural complexity of neutral atoms using the non-relativistic 
Hartree-Fock (HF) wave functions 49|| for atoms with atomic number Z=l-54, have been 
reported using a variety of information measures. A similar evaluation of a complexity mea- 
sure for neutral atoms with Z=l-103 was recently carried out in terms of the Dirac-Fock 
wavefunction, choosing the exponential of the Shannon position entropy as the measure of 
information, H = exp(S'r), and the average position density as the measure of disequilib- 
rium, D =< p > jsl. While the measure of information exhibits a strong shell effect but 
insignificant relativistic effect, the measure of disequilibrium was found to be a monotoni- 
cally increasing function of Z exhibiting a strong relativistic effect. It is remarkable that the 
ratio of the relativistic to the nonrelativistic measures of disequilibrium (average position 
densities) obtained for neutral many-electron atoms is in almost quantitative agreement with 
the corresponding ratio obtained in terms of the average position densities evaluated above 
for the single electron ions. These ratios are presented in Fig. 4. This observation must be 
a manifestation that the relativistic effect on the measure of disequilibrium is dominated by 
the effect on the innermost orbital. 



VII. SCALE INVARIANT ENTROPIES 

A. Residual position and momentum entropies 

The KuUback-Leibler relative/residual information measure of a given probability density 
is defined with respect to a prior density and it determines the extra information contained 
in the given density relative to the prior. Such a residual entropy for the relativistic density 
over the nonrelativistic density as the prior can be defined in both the position and the 
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momentum space. In position space this can be done analytically, yielding 



ATrr'pRir) log J dr (35) 

^ 0.197467(aZ)^ + 0.150105(aZ)^ + 0.115104(aZ)^ + • • • 

The Taylor series for ^5^^ was obtained by evaluating its first six derivatives with respect 
to Z, using maple. The values of S?^^^ for Z < 25 were calculated using the three-term 
Taylor series, since evaluating the analytic expression involves cancellation errors. We note 
that for Z = 25 the analytic expression and the three-term expansion practically coincide. 

The extreme relativistic value is obtained from eq. (|35il by evaluating its limit as 7 0. 
It is found that Sr^^^ = log(2) + 2C ^ 1.847579. 

The residual momentum entropy Sp^^^ = Airp'^Ilji^p) log ^ yinr{p) ) evalu- 
ated numerically. Differentiating the integrand four times with respect to Z and inte- 
grating numerically we obtained the leading term in the power series expansion Sp^^^ ~ 
4 



0.572467(aZ)^ + ■ 

s 

(aZ) 



gR/NR 

The value of / „^4 at Z = - was obtained using the extreme relativistic momentum 



density, eq. fl20l) . 

Since Sr ^ and Sp ^ are pure (dimensionless) quantities, they do not measure position 
or momentum widths or uncertainties. They do measure the (somewhat slow) onset of 
relativistic effects upon increase of the nuclear charge. 



B. Average Measures of relative distance 



The average measure of relative distance 28|, |50| of the position densities is the sum of 
the two relative entropies 



SR/NR = / Anr'pnir) log . , , 
Jo \PNR{r) 



(_Pr{ 



dr 



and 



S^^l^ = / Anr'pr^Rir) log ( ^^^^ ) dr 

Jo \ pRKn 
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It can be written in the form 



The measure of relative distance of the momentum densities is defined in an analogous 



manner 

gR/}^R was evaluated above, cf. eq. (l35l) . Sr^^^ can be evaluated in a similar way, 
yielding 

Sr^^ = log (^^^^^) + (1 - 7) (3 - 2C) . 

Hence, 

= (1 - 7) ( 3 - 2C - 2^(27) ^ 



7. 

The Taylor series for 5*^ can be obtained analytically. The first three terms are given by 
S.«(^-5)(.zr.(gH-C(3)-g(.Z)« 

V90 96 4^^ ^ 64 y ^ ^ 
^ 0.394934(aZ)^ + 0.27452(aZ)'^ + 0.20103(aZ)^ + ■ ■ ■ 

Sp was evaluated numerically. 

The residual (relativistic vs. nonrelativistic) position and momentum entropies, and the 
average measures of the distances of the corresponding position and momentum distribu- 
tions, are presented in Table 1, all normalized via division by (aZ)^. We note that is 
a monotonic function of Z, but -t^Stt is not. 



VIII. CONCLUSIONS 



The characterization of inherent quantum mechanical uncertainties has become a rich 
field of study with direct relevance to emerging technologies. In the present article we 
examine the application of widely used information measures to the ground state of the 
relativistic hydrogen-like atoms, clearly bringing out the dependence on Z due to the rela- 
tivistic effects. Further, we point out and illustrate the well-established but largely ignored 
difficulties associated with the most common quantum mechanical formulation of the un- 
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certainty principle, that arise as a consequence of the fact that the radial momentum is 
not self-adjoint. Several information measures exhibit singularities at particular nuclear 
charges, notably ^ = 2^ ~ 118.68 and Z = ^ 132.68, whose significance remains to 
be elucidated. In the coordinate representation all the information measures considered 
allowed analytic evaluation of the integrals involved. This has not been the case for the cor- 
responding momentum space quantities. What we find particularly puzzling in this context 
is the fact that the closed analytic expression for the position-space expectation value of the 
Laplacian agrees, as expected, with the numerically evaluated average over p^, in momentum 
space, and still we failed to evaluate the latter analytically. These, and many other issues 
such as uncertainty and information measures for excited states as well as for many-electron 
atoms, suggest that the study of information measures for relativistic systems is a widely 
open field. 
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TABLE I: Residual entropies and average measures of relative distance. 
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gR/NR 


gR/NR 


(aZ)4 
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3.000000 


5.485774 



34 



Momentum 




Z 



FIG. 1: Ratios between the relativistic and nonrelativistic Shannon and Fisher lengths {Ps^Pf) 
and impeti (^5, /ii?), and corresponding ratios for root mean square of position (/?) and momentum 
(/u) as functions of Z . 
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FIG. 2: Ratios between the relativistic and nonrelativistic Renyi lengths (Pr) and impeti {ij,r) as 
functions of Z. 
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FIG. 3: Renyi sum vs s for the non-relativistic H atom and relativistic H-like atoms. 
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FIG. 4: Ratio of relativistic to non-relativistic estimates of the linear entropy for neutral atoms 
and H-like atoms. 
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